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Let  US  cHcck  the  functionals  (39)  and  (40)  by  studying  their  special  cases  for  the  isoiropit  guide,  i  r  ,  tor  f  •  f  0  Wnung  L 
(fit  -  i>^)%and  =  M ,  wt  have 
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This  expression  can  be  seen  to  coincide  with  K^at  given  in  {lOj 
(note  that  the  factor  1  / in  the  (21)  of  the  refe>Mccd  paper  is 
erroneously  missing). 

III.  Discussion 

The  functional  (39)  is  exact  and  presents  an  obvious  cxtcnsibo 
of  previously  known  functionals  from  isotropic  to  bi-isoiropic  me-^ 
dia.  It  can  be  applied  to  obtaining  approximate  mode  solutions  for 
open  bi-isotropic  waveguides.  By  approximating  the  longitudinal 
held  functions,  the  dispersion  function  /3(u)  can  be  obtained  point 
by  point.  In  practice,  this  can  be  performed  using  the  following 
scheme: 

1 .  Choose  a  value  for  the  phase  velocity  parameter  t(p. 

2.  Find  suitable  approximating  functions  for  the  longitudinal 

fields  e(  p)  and  h{  p),  i.e.,  the  matrix/(  p),  with  free  param¬ 
eters.  Insight  on  the  field  distribution  of  the  problem  will 
help  in  allowing  use  of  just  a  few  parameters;  otherwise,  a  , 
massive  computation  scheme  with  a  gteat  number  of  patany/ 
eters  is  needed.  / 

3.  Optimize  values  of  these  parameters  so  that  the  func^nal 

HVp,  f)  obtains  the  stationary  value;  i.e.,  its  differ^iation 
with  respect  to  all  these  parameters  is  zero.  (In  ca^of  large 
number  of  parameters  this  requires  use  of  some  wtimizat'on 
procedure.)  / 

4.  The  corresponding  parameter  values  inserte^in  the  longitu¬ 
dinal  field  expressions  give  closest  appmimations  to  the 
fields  and  the  functional  value  approxim^es  the  value  of 

5.  Now  it  is  ea^to  determine  a  point  fln  the  dispersion  dia¬ 
gram:  oi  =  vu)^,  0  —  wVp.  / 

0.  The  transverse  field  functions  corrKponding  to  this  point  are 
obtained  from  (23).  / 

7.  To  find  another  point,  start  vvjui  another  value  of  Vp. 

Thus,  the  procedure  works  best ^hen  some  a  priori  knowledge  of 
the  longitudinal  fields  exists,ywhich  helps  in  finding  suitable  ap¬ 
proximating  functions  withal  too  many  optimizable  parameters. 
Obviously,  the  method  is^pecially  attractive  for  finding  the  low- 
est-order  modes  with  le^  spatial  variation.  To  find  the  knowledge 
required,  it  appears  n^ssary  to  work  through  some  examples  with 
brute-force  technique.  This  is,  however,  outside  the  scope  of  the 
present  theoreticaj/study 
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Analysis  of  Bilateral  Fin-Lines  on  Ani.sotropic 
Substrates 

Thinh  Quoc  Ho  and  Benjamin  Beker 


Abslracl—A  full-wave  analysis  of  the  bilateral  fin-line  on  anisotropic  j 

substrates  is  presented.  The  supporting  medium  is  characteriied  si¬ 
multaneously  by  both  nondiagonal  second  rank  (e)  and  |  ul  tensors.  The 
dyadic  Green's  function  is  formed  rigorously  in  the  discrete  Fourier 
transformed  domain  and  is  used  to  study  the  propagation  character¬ 
istics  of  the  fin-line.  The  Green’s  function  elements  are  given  explicitly 
in  their  closed  forms  along  with  the  verification  of  the  theory.  New  data 
describing  the  dispersion  properties  as  functions  of  the  Cfwrdinale  mis¬ 
alignment  are  also  generated  for  several  substrate  materials. 

I.  Introduction 

IQ  CtU 

Although  the  theories  of  transmission  lines  on  anisotropic  stnic  fnjP 
lures  are  well  documented,  the  major  effort  thu,s  far  has  been  di  ** 


Manuscript  received  June  3.  1991 .  revised  October  4.  ]9>) 

The  authors  arc  with  the  Dcpartmcni  of  l-iccinca!  and  <  omputcr  F.ngi 
neenng.  University  of  South  Carolina.  Columbia,  SC 


IEEE  Log  Number  910477^ 


0018  9480/92503  00  0  1992  IEEE 


406 


!H  t  TRANSACriONS  ON  MtCKOWAVfc  fHKiKV  AM)  IK  HNlyul  S  '.ill  4H  Nl)  1'  tlHRl  AK)  i«i 


reeled  toward  structures  whose  subsitai-s  are  eharacien/cd  by  a 
dtagonalaed  (tj  tensor  only  A  review  ot  the  ircalmenis  to  those 
problems  outlining  the  procedure  for  the  study  of  the  propagation 
characteristics  of  various  transmission  lines  may  be  found  in  11] 
|3)  In  practice,  however,  most  of  the  abtjve  analyses  arc  somewhat 
limited  in  scope  due  to  the  absence  of  the  off  diagonal  elements  in 
the  permittivity  tensor,  which  may  be  used  to  represent  the  mis 
alignment  between  the  matcnal  coordinate  system  and  that  of  the 
waveguide,  as  indicated  in  |4)  In  addinon,  the  respsmse  to  the 
magnetic  fields  for  the  aforementioned  cases  is  normally  assumed 
to  be  isotropic,  that  is,  when  |mI  is  specified  as  a  ^ero-rank  tensor 
(scalar).  For  problems  involving  transmtsston  lines  on  general  an¬ 
isotropic  substrates,  that  is,  when  the  materials  are  characterized 
by  both  permittivity  and  permeability  tensors,  the  solution  to  Max 
well’s  equations  required  for  the  study  of  their  propagation  char 
acteristics,  can  no  longer  be  obtained  in  a  straight  forward  manner 
Up  until  now,  there  are  only  a  few  publications  I5|-(S|  which  have 
reported  some  research  efforts  dealing  with  transmission  lines 
whose  substrate  material  is  characterized  simultaneously  by  both 
[e]  and  |m)  tensors.  The  fact  that  the  shielded  transmission  lines  on 
a  general  anisotropic  media  have  not  yet  received  full  attention,  as 
did  their  counterparts  printed  on  substrates  characterized  by  (<!  ten¬ 
sor  alone,  provides  the  motivation  for  this  work 
The  purpose  of  this  paper  is  to  present  the  analysis  of  bilateral 
fin-lines  on  general  anisotropic  substrates,  in  this  study,  unlike 
some  others  conducted  in  the  past,  both  [e]  and  (p)  tensors  are 
nondiagonal.  As  a  result,  the  Green’s  function  denved  herein  may 
be  used  to  examine  the  dispersion  properties  of  the  structure  on 
both  dielectric  and  magnetic  substrates  The  bilateral  fin-line  is  se¬ 
lected  mainly  due  to  its  popularity  in  applications  at  millimeter- 
wave  frequencies.  Also,  this  type  of  a  transmission  line  docs  not 
support  a  TEM-ltke  mode.  It  has  a  cut  off  at  lower  frequencies,  and 
as  such  is  fundamentally  different  from  other  transmission  tines  At 
the  present  time,  there  is  no  comprehensive  treatment  of  this  struc¬ 
ture  that  is  printed  on  the  type  of  anisotropic  materials  considered 
in  this  study.  Consequently,  the  aim  of  this  work  is  to  present  the 
analytical  tools  which  can  be  used  in  the  analysis  of  bilateral  fin- 
lines  whose  substrates  can  exhibit  a  wide  variety  of  anisotropic 
properties. 


II.  Theory 

The  cross  .section  of  the  bilateral  fin-lmc  is  show  n  in  Fig  1  along 
with  the  coordinate  system  used  in  the  analysis  The  millimeter 
wave  circuit  using  the  fin-line  is  usually  excited  by  a  dominant 
TEio  mode  of  the  rectangular  waveguide  having  dimcn.sions  <i  and 
b.  The  direction  of  propagation  of  the  resulting  wave  is  along  the 
Z-axis.  In  this  geometry,  the  gap  that  .separates  the  metal  fins  is 
denoted  as  W  The  supporting  substrate  oniu  which  the  tins  are 
mounted  is  suspended  along  the  fj-plane  ol  the  waveguide,  and  is 
specified  by 
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that  may  exist  within  the  planar  anisotropic  region,  the  two  curl 
equations  arc  employed  It  can  be  .show  n  that  the  vector  wave  equa 
lions  for  all  components  of  the  electric  and  magnetic  fields  iiiav  l>f 
written  compactly  in  the  form  of 

VX(Ip)  ‘  VXff)  H-jil  F:  0  1 2a I 

VX([<|  '  rx/f)  i'liij  •  //  0  (2b) 

where  is  the  free-spacc  wave  number 
To  obtain  the  field  solutions  within  the  anisotropic  medium, 
cither  the  wave  equation  for  £  or  the  wave  equation  for  H  may  be 
used.  It  is  more  convenient  to  work  with  the  wave  equation  for  the 
electric  field  mainly  because  the  problem  involves  several  ciccinc 
ground  planes  which  belong  to  the  housing  Consequently,  the  vet 
lor  wave  equation  for  the  electric  field  may  be  reduced  lo  ihice 
scalar  equations  which  can  be  written  in  term.s  of  the  cletiiu  field 
components  alone.  The  priKcss  of  separating  (uncouplingi  these 
equations  to  get  a  set  of  independent  cqualmns  for  individu.it  i  om 
ponenis  of  the  £■  field  begins  with  then  transformation  to  the  1  ou 
ncr  domain 

The  spectral  representation  ol  any  field  coirifvmeni  is  di  fiiied  \i,i 
the  following  transform 
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Z,  =  a&liyy/liy,  -  „  j  t  ^  -  *»<..)Z„.  (5k) 


These  coefficients  are  functions  of  the  medium  parameters,  trans¬ 
form  variable  or,  t„.  and  tl  c  propagation  constant  /S  which  is  de¬ 
fined  along  the  z  direction. 

Decoupling  of  the  above  set  is  possible  by  the  substitution  method 
which  after  some  manipulations  leads  to  yet  another,  but  an  inde¬ 
pendent  pair  of  fourth-onJer  differential  equations  for  these  com¬ 
ponents  of  the  electric  field,  i.e.. 


a"!?,  a^f,  a^£,  a£, 

TT  +  ^  T-t  +  B  -rr  +  C-—  +  D£,  =  0 
die*  ax’  dx^  dx  ’ 


(6a) 


d*£.  d'£,  d^£.  d£. 

where  the  transformed  constants  /4,  fl,  C,  and  D  are  defined  below: 

TjZj  +  T,(n  -  K3Z,)/K, 


4  =  -  ,n  -  w,,/r.  4 

B  =  Y,  -  y,Z,  -  KifT,  -  T,z,)/)',  -E  (T,  -E  YA/Y,) 

Y,Z2  e  Y,(Y4  -  Y,z,)/Y, 

Y,  -  Yy{2,  -  ZtYy/Y,) 


(7a) 


(7b) 


C  =  -r,Z.  -  Y2(Y4  -  T,Z,)/K, 


E 


iY,  y,z,  E  y,z,yjy,) 


y,7^  e  y.(y,  -  y^z,)/y, 
n  -  r,(Z,  -  Zjy,/^,) 


D  -  -  y,(Z4 


V  7  y ,  ZtZ;  E  n(K.  -  nz,)/y, 
n  -  y,(Z,  -  z,y,/Y4) 


(7c) 

(7d) 


and  where  the  coefficients  E,  F,  G,  and  H  arc  also  functions  of  K,, 
^2-  Y4,  Z|.  Zj.  Z),  and  Z^.  Within  a  shielded  housing,  the  so¬ 

lutions  to  (6)  are  standing  waves.  In  general,  these  solutions  arc 
written  in  terms  of  sinusoidal  and  cosinusoidal  functions.  Due  to 
the  symmetry  of  the  structure,  the  boundary  conditions  require  the 
tangential  magnetic  fields  to  be  zero  at  the  magnetic  wall  (jr  '  0), 
reducing  the  general  solutions  to  having  cosine  functions  only 
Subsequently,  when  the  two  tangential  electric  field  components 
have  been  determined,  the  remaining  field  components  such  as 
£,.  /?. .  W. .  and  H.  may  be  expressed  in  terms  of  £,  and  £.  through 
the  Maxwell's  curl  equations.  On  the  other  hand,  in  the  isotropic 


region,  the  field  cotnponents  iiuv  tw  found  ()y  using  ifit-  scaiji  po 
icnttal  functions  I  he  procedure  (or  tindirig  them  is  aiie.>dv  jsji! 
able  in  (9).  and  therefore,  will  not  t>e  re|)eatcd  here 

Once  all  (he  field  conqKinciits  m  txjth  regions  ate  know!!,  the 
Green's  function  of  the  struciure  may  (>c  obtained  by  enlort  lag  the 
appropriate  boundary  conditions  at  the  air-subsiratc  interlace 


■  £,: 

(ha 

7.'; 

shb 

H,: 

iH-. 

tl,  t'l,,  -■  J, 


(Wi 


where  J,  and  J.  are  the  l-'ounct  transforms  of  unknown  correm 
densities  J,  and  J.  an  the  fins  at  t  -  h, 

After  some  mathematical  manipulations,  (he  dyadic  Green's 
function  in  the  Founcr  iransfonned  domain  is  found  to  relate  the 
fields  across  the  slot  to  the  currents  on  the  fins  through  the  follow¬ 
ing  malnx  equation 


[C] 


with  the  expressions  for  G„.  (J.,,.  given  m  the  Appendix 

Notice  that  £.  and  £^  in  the  above  equations  are  finite  across  the 
slot  at  a  =  A,  for  0  <  j  y\  <  w/2,  but  they  are  zero  on  the  metal 
fins  In  order  to  find  the  numerical  solution,  the  apenure  fields  are 
expanded  using  the  known  basis  functions  19), 

u 

£,(n)  =  Z  C.£»(cr(n))  (lOa) 

^  •  -  I 

H 

£.(n)  =  21  D.,£.Jnin))  (10b) 

•  -  1 

with  and  D„  being  the  unknown  expansion  coefficients  After 
the  appropriate  substitutions  between  (9)  and  (10)  along  with  the 
application  of  the  Parscval’s  theorem  in  the  discrete  Fourier  do¬ 
main.  the  system  of  matrix  equations  which  is  used  to  extract  the 
propagation  constant  8  is  formed 

M  V 

^  E  V  AU  ’D„  0  ,  .--  I.  2,  7.  -V 

•f  I  «»  X.  } 

(Hat 

M  ,V 

2-  a;„”C,  e  E  a'.  D.  -  0  I  '  I.  2.  7.  .V 

•  I  m  ■"  I 


(lib) 

with  Its  elements  ei  .^cn  by 

K,Jp.q)=  E,^lfi}C  In.  It)  F,Jnj  flit) 

"  0 


where  :>  or  q  can  be  either  v  or  ;  The  roots  which  correspond  to 
the  Tiropagation  constants  arc  found  by  setting  the  determinant  ol 
the  coefficient  matrix  equal  to  zero,  and  solving  the  resulting  eqiia 
lion  (11)  using  the  technique  described  in  (9) 


III  Rfstii  IS 

'll)  valKiaic  the  newly  derived  expressions  of  the  Grceit'-.  tunc 
lion  elements,  numerical  results  for  the  effectisc  diclccttK  cunstani 
(lt  /l;i  far  two  special  cases  are  (ompiiicd  and  lomp.ircd  lo 
the  existing  data  In  the  first  case,  ific  supporting  mednini  i'-  i,iken 
lobe  isotropic  wilhi,,  i,,  <  i  'S.  u,,  n,,  u  10 
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and  <„  =  =  M.,  =  A*,.  =  0.  The  chosen  housing  is  a  slandard 

WR-28  waveguide,  with  the  remaining  dimensions  given  by  h,  - 
0.0625  mm,  =  3  4935  mm.  and  =  0  5  mm  For  tlie  second 
case,  all  physical  dimensions  of  the  structure  are  the  same,  but  the 
substrate  is  a  dielecincally  anisotropic  medium  with  «„  -  t  ..  - 
3.0,  f„  -  3.5,  M,,  =  Mv,  =  A*;:  =  i  0-  and  t,.  =  -  ji.,  =  a*,. 

=  0.  Fig.  2  shows  the  effective  Jielecine  constant  computed  by 
this  melhtxl  along  with  the  data  produced  from  iZ],  19].  The  com 
panson  of  the  results  indicates  that,  in  general,  the  agreement  is 
very  good.  For  the  isotropic  substrate,  the  1,^  matches  very  well, 
including  frequencies  below  the  designated  bandwidth  (26.5  GHr. 
40.0  GHz).  However,  for  the  uniaxial  material,  a  minor  discrep 
ancy  is  observed  at  frequencies  near  the  cut  off  point,  but  other¬ 
wise,  everywhere  else,  the  agreement  is  good  This  difference  in 
the  results  near  the  cut  off  can  probably  be  attributed  to  the  nu¬ 
merical  convergence  of  the  data  obtained  by  the  two  methods.  In 
all  of  the  calculations  presented  in  this  paper,  the  values  of  N  =  S. 
A/  =  5,  and  n  =  350  terms  were  used  to  ensure  convergence. 

Now  that  the  Green’s  function  has  been  validated,  the  effects  of 
the  nondiagonal  elements  in  the  [e]  and  (/i|  tensors  on  the  disper¬ 
sive  properties  of  bilateral  fin-lines  can  be  examined.  First,  the 
permeab.''ty  tensor  of  the  substrate  is  made  isotropic  by  allowing 
Mxj  =  M,y  —  Uzz  =  10  Miy  and  a*,,  =  0.  The  effective  dielectric 
constant  of  the  fin-line  is  then  studied  as  a  function  of  the  rotation 
angle  $  which  in  practice  may  be  used  to  represent  the  misalign¬ 
ment  between  the  coordinate  system  of  the  waveguide  and  that  of 
the  substrate  material.  This  leads  to  the  following  definitions  for 
the  tensor  elements  of  the  permittivity  appearing  in  (la):  <„  =  <2 
sin^  (tf)  +  «i  cos^  (ff),  e„  -  <2  (^)  +  fi  sin^  (&).  and 

-  (<J  “  <i)  S'"  (®)  cos  (6).  where  t,.  «2,  and  <3  are  the 
principal  values  of  the  tensor.  Numerical  results  in  Fig.  3  show 
the  response  of  <,»  when  the  fins  are  primed  on  the  PTFE  or  glass 
cloths,  both  of  which  are  dielectrically  biaxial  substrates.  The  ef¬ 
fective  dielectric  constant  of  the  fin-line  is  computed  at  frequencies 
of  26.0  GHz.  33.0  GHz,  and  40,0  GHz  and  is  plotted  versus  the 
rotation  angle  S.  As  can  be  seen  from  this  figure  (^(r  becomes  larger 
with  increasing  frequency,  however,  for  both  of  the  cloth  materials 
it  is  decreasing,  but  only  slightly,  as  the  angle  &  varies  from  0  to 
90  degrees.  This  behavior  is  totally  opposite  to  that  of  the  micro- 
strip  line  wherein  it  is  found  e,iT  increa.ses  with  values  of  the  rota¬ 
tion  angle  changing  from  0°  to  90°,  as  illustrated  in  (lOJ. 

Next,  the  response  of  t,,f  belonging  to  the  fin-line  structure  which 
uses  a  dielectrically  uniaxial  substrate  is  examined  The  two  chosen 
substrate  materials  are  boron  nitride  and  sapphire  with  the  physical 
dimensions  of  the  housing  being  the  same  a.s  tho.se  u,scd  above.  A 
very  interesting  behavior  can  be  observed  from  the  numerical  re¬ 
sults  which  are  plotted  in  Fig.  4  In  the  case  of  a  fin-line  printed 
on  sapphire,  the  value  of  the  effective  dielectric  constant  increases 
with  the  rotation  angle,  while  for  the  boron  nitride  substrate  the 
opposite  is  true. 

Finally,  to  demonstrate  the  flexibility  of  the  theory  presented  in 
this  paper,  the  propagation  properties  of  bilaicral  fin-lincs  using 
more  general  anisotropic  media  arc  computed  1  he  medium  perme¬ 
ability  is  now  characterized  by  a  second  rank  tensor  whose  cic 
ments  arc  given  by  /x,,  =  ft,  sin’  (y:)  t  p,  cos'  (y?).  rt,,  =  ji;  cos’ 
(vi)  +  Ml  (v>)-  M:.-  =  M}-  atlAl  M,,  M,,  '  <  M.’  '  Alt)  sin  (y?)  cos 

(v>).  In  this  case,  variables  <p.  Mi-  A'.’  and  /i,  are  the  anatoges  of  0. 
fi,  (j,  and  fj,  respectively.  Fig  5  illustrates  how  the  normalized 
wavelength  X,  reacts  to  the  addition  ol  the  magnetic  anisotropy 
The  substrate  material  parameters  used  in  computations  were  r,  - 
8  7,  <2  =  9.6,  (\  =  I  I  4.  Pi  -  I  0.  p  .  1  h.  and  p,  -  I  8  The 

calculated  data  indicates  that  the  [iiesense  ol  the  |p|  tensor  ctfec 
lively  shortens  the  normalized  wavelength  ol  the  lin  line  consul 
crably  This  effect  can  be  seen  lot  hoiti  diagon.il  ip-  D”.  90°)  as 
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Fig  2  Bffcciivc  diclcciric  consUfif  of  (h<r  biiaterrai  hn  hnc  a  •-  1  sl2 
mm.  b  -  }  556  mm.  A,  =-*  0  (>625  mm.  h;  ^  49,^5  mm.  -  0  5  mm  on 
ivijnopic  substrate  3  75  and  ^  /«,  -  I  0.  and 

dielecincally  unitmal  substrate  wjth<,  -  >  0.  «  3  5.  and  tt.  - 

-/*,-)  0 
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Fig  3  FfTcctjvc  diclcctnc  constant  versuv  B  of  the  bilateral  hn  line  on 
PTFF  and  glass  cloths  with  a  '■  ?  112  mm  h  3  556  mm.  h,  0  (>62'' 
mm.  A;  -  3  49.35  mm,  h-  0  s  mm 
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Fig  4  F.Hcclivc  diclcctfu  (.<M3\f,in[  vrrMf's  $  id  the  hiljtcral  tin  line  on 
boron  nitndc  and  sapphire  kkjiti  d  ^  !  1.'  mm,  b  ^  ^56  mm. 
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well  as  thHi  diagonal  (v"  •1"'  >  pofnicahilils  tensor,  sshrtcin.  lor 

Ihis  set.  all  results  were  i>hfaiiu’ii  .ii  the  (.cnior  (u'ljucocv  of  (he 
At;  band 
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Fig.  5.  Nonnalized  wavelength  versus  9  as  a  function  of  ^  of  the  bilateral 
fin-line  on  general  anisotropic  substrate  with  a  =  7.1 12  mm.  b  =  3.5S6 
mm,  ft,  =  0.0625  mm,  Aj  =  3.4935  mm,  w  =  0. 15  mm,  r,  =  8.7,  «,  = 
9.6,  *j  =  11.4.  Ml  =  l.O.  mj  =  1 .6.  and  jtj  =  1 .8,  and /  =  33.0  GHi.  (a) 
(e]  only,  (b)  [<)  and  [/»]»>  =  0*  (c)  («1  and  (Mis’  =  45*.  (d)  (<)  and  ImIv 
=  90“. 


IV.  Conclusion 

A  fiill-wavc  analysis  based  on  the  spccfral-domain  technique  of 
bilateral  fin-lines  on  anisotropic  substrates  was  presented.  The  for¬ 
mulation  of  the  Green's  function  was  carried  out  in  the  discrete 
Fourier  transformed  domain.  Both  the  magnetic  and  dielectric  an¬ 
isotropy  effects  on  bilateral  fin-lines  were  examined  to  illustrate  the 
usefulness  of  the  newly  derived  Green's  function  for  this  structure. 
The  dispersion  properties  were  computed  for  several  substrate  ma¬ 
terials  as  a  function  of  the  misalignment,  and  it  was  found  that  in 
some  cases  they  can  be  greatly  affected  by  the  rotation  of  the  co¬ 
ordinate  axes.  These  generalized  expressions  may  also  be  used  to 
study  the  propagation  characteristics  of  other  transmission  line 
structures,  besides  the  bilateral  fin-line,  when  they  are  printed  on 
dielectric,  magnetic,  as  well  as  mixed  anisotropic  media. 


Appendix 

The  dyadic  Green’s  function  elements  appearing  in  (9)  arc  de¬ 


fined  below 

C„(a,  ff)  =  6,/*^  +  {//;+*  -  H:*~)/62  (Ala) 

G,^(a.  p)  =  63/t^  +  +  e;h:9-}/52  (Aib) 

&)  =  6^/*^  +  -F  (Ale) 

G,,(0f,  0)  =  -F  {£;;,v;'f‘  -  e:h,  '^  ]/5,  (Aid) 

♦  -  =  tan  (yjfij),  (Ale) 

♦  *  =  tan  (y*/i|),  (Alf) 

♦  =tan(7.)r,)  (Alg) 

with 

=  (*i  -  0^)/.y7  (A2a) 

6,  =  (E*  -  £,:)/x,-  (A2b) 

6,  =  a0/y2  (A2i) 

l>4  -  (Ey  -  f,' )nvv  (A2d) 

6,  =  {kl  -  a')/7j  (A2i.-) 


£■;=  F  zlY.ta: 

t  /.yZlYs)] 

+  {-2A'  -  (Zinz.  -  zinzi) 

iZl  -  l\7.sYi  F  ZsZlYi)] 

(A3a) 

Er  =  b  +  zlY,)a: 

ZU  Y,  F  Za/iTj)! 

+  i-ZsZ,'  --  a\YA  -  2.\Ys7.s) 

al  -  ZlZsYi  F  Z^ZlY,}) 

(A  3b) 

=  {I/l-jiSy.d  - 

'  (yI  ~  F  a‘iir,/i‘u 

"  E*  a0n„/is^ 

+  Fjr  Ci0tS  fy/ fSy^(  ..f} 

(A3c) 

He  =  {l/l-J0Y-(l  -  »»')!} 

'  {>-  - 

a: /Mutf.  -  Ej  a0is„/ft,y 
+  Ey  a0ii„l„/2i^(„) 

(A3d) 

H:  =  -  jcH„  +  yyE: 

(A.3e) 

H;  =  ~jaHg  +  y.Ey 

(A3f) 

H*  =  Pr,MuPa'(-Ya  +  j0HJ 

(A3g) 

Hr  =  PrrMuPa'(-7-  +  j0Hg) 

(A3h) 

where  72  is  the  transverse  wavenumber  in  the  isofropic  region.  Pa¬ 
rameters  7^  are  the  corresponding  wavenumbers  in  the  anisotropic 
region,  and  they  are  obtained  by  solving  the  fourth-order  charac¬ 
teristic  equation  using  a  technique  similar  to  the  one  presented  m 
IIOJ. 
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